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Abstract: Nuclear security usually requires the simultaneous detection of neutrons and 
gamma rays. With the development of crystalline materials in recent years, Cs2LiLaBre 
(CLLB) dual-readout detectors have attracted extensive attention from researchers, 
where real-time neutron/gamma pulse discrimination is the critical factor among 
detector performance parameters. This study investigated the discrimination 
performance of the charge comparison, amplitude comparison, time comparison, and 
pulse gradient methods and the effects of a Sallen—Key filter on their performance. 
Experimental results show that the figure of merit (FOM) of all four methods is 
improved by proper filtering. Among them, the charge comparison method exhibits 
excellent noise resistance; moreover, it is the most suitable method of real-time 
discrimination for CLLB detectors. However, its discrimination performance depends 
on the parameters t,, tm, and te. When t, corresponds to the moment at which the 
pulse is at 10% of its peak value, tẹ requires a delay of only 640 to 740 ns compared to 
t,, at which time the potentially optimal FOM of the charge comparison method at 
3.1-3.3 MeV is greater than 1.46. The FOM obtained using the t,, value calculated by 
a proposed maximized discrimination difference model(MDDM) and the potentially 
optimal FOM differ by less than 3.9%, indicating that the model can provide good 
guidance for parameter selection in the charge comparison method. 
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1 Introduction 


Neutron and gamma detection is an essential technique for many nuclear security 
applications, such as isotope identification '7!, radiation monitoring P451, and the 
search for fissile material [67], The detection of neutrons and gamma rays usually 
requires a combination of both types of radiation detectors, which must be portable and 
compact in many detection scenarios. In recent years, a new Cs2LiLaBre (CLLB) 
detector in the elpasolite scintillator family has been developed '!. It has the 
advantages of small size, high light output, and excellent energy resolution; it enables 
the detection of two types of radiation using pulse shape discrimination technology, 
making handheld instruments possible !"'7!, Since the 1960s, researchers have found 
that scintillators produce different pulse shapes when they interact with different 
particles |!3], and many methods of discriminating the waveform characteristics have 
been developed. These methods are divided into three categories: time-domain, 
frequency-domain, and intelligent methods. Time-domain methods include charge 
comparison !'4! (CC), time comparison (TC), amplitude comparison (AC), and pulse 
gradient (PG) methods !!51, Frequency-domain methods include frequency gradient !"°, 
wavelet transform |'7!, and fractal spectrum !'*! methods. Intelligent methods include 
multilayer perceptron l°, support vector machine ?°"), and deep learning network 7! 
approaches. Time-domain methods are simple to calculate but usually have a low 
figure of merit (FOM). Frequency-domain methods involve many Fourier calculations 
that are complex and difficult to run in real-time; however, they are less sensitive to 
noise than time-domain methods. Intelligent methods usually have the highest FOM 
but require many floating-point matrix operations and thus a complicated 
computational process. In addition, the discrimination performance of intelligent 
methods usually depends on the training set size and number of neural network 
parameters. It is almost impossible to complete the computations required by 
frequency-domain and intelligent methods using limited computational and storage 
resources. In addition, as the sample points of the waveform increase, the 
computational complexity of frequency-domain and intelligent methods can rarely 
maintain linear growth. Therefore, time-domain methods are still preferred for real- 
time analysis. Consequently, it is necessary to study suitable discrimination methods 
and preprocessing techniques to identify methods with a good FOM and some noise 
immunity that are suitable for real-time analysis using modest computational resources. 


Zuo et al. [I3] explored the discrimination effect of several time-domain methods on 


plastic scintillator detectors and found that filtering could enhance the discrimination 
performance of these methods. However, the applicability of this finding to CLLB 
detectors is unknown. In addition, time-domain methods are all sensitive to their 
parameters. To optimize the discrimination performance of these methods, one must 
continuously search for the appropriate parameters by trial and error, which is very 
laborious. 

To solve these problems, we examine the discrimination effects of the CC, TC, 
AC, and PG time-domain methods on CLLB detectors and the effects of preprocessing 
by a Sallen—Key (SK) filter on the methods in real time. A new SK recursive 
expression based on the second-order Runge—Kutta method is proposed to reduce the 
error and applied in an experiment. The experimental results show that appropriate 
filtering can enhance the discrimination performance of these time-domain methods, 
where the CC method is the most suitable real-time analysis method for CLLB 
detectors. However, the performance of the CC method is affected by the parameter 
settings. We discuss the effects of the CC parameters on its discrimination performance 
and propose a maximized discrimination difference model (MDDM) to guide 
parameter selection in the CC method. This tool can be used to model the parameter 
settings of other time-domain methods and guide parameter selection. 


The remainder of the paper is organized as follows. The experimental setup for 
the work is described in Sect. 2. The basic principles of the numerical methods used 
are described in Sect. 3. The discrimination effects of the methods and the effects of 
filtering and the MDDM on their performance are discussed in Sect. 4. Conclusions are 


presented in Sect. 5. 
2 Neutron/gamma pulse discrimination experiment 
2.1 Detector and radioactive source 


A CLLB crystal detector from Saint-Gobain '*! was used; the scintillator crystal 
has a diameter of 2 in. and is sealed inside an aluminum housing. It was coupled 
directly to a photomultiplier tube (PMT) (Hamamatsu R6231-100), as shown in Fig. 1. 
An ORTEC 556 high-voltage power supply was used during the experiment to supply 
a voltage of 800 V to the detector. A '°’7Cs gamma source was used as the energy 


calibration source for the detector, and a *’Cf source was selected to provide the 


neutron/gamma hybrid radiation field. 
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Fig. 1 Neutron/gamma signal discrimination system based on CLLB detector 


2.2 Data acquisition system 


A data acquisition system developed by our research group was used for the 
experiment; it includes a PXIe chassis with an embedded controller and an electronics 
card. The signal input of the main electronics card was directly connected to the CLLB 
detector. The main electronics card consists mainly of analog signal conditioning 
components, high-speed analog-to-digital converters (ADCs) (14 bits and 500 Msps), a 
clock jitter cleaner, and a field-programmable gate array (FPGA) (type 
xc7k325tffg900). As shown in Fig. 4, the rise time of the electrical signal of the CLLB 
detector is t, = 30 ns. If the signal is processed directly by a 500 Msps ADC, there are 
no fewer than 15 sampling points on the rising edge. A signal with a rising edge t, = 
30 ns corresponds to the frequency spectrum of 16.7 MHz, which is much lower than 
the sampling rate of the waveform digitization module for 500 Msps ADCs. The 
hardware of the main electronics card is almost identical to that of the main electronics 
card for neutron flux monitoring in the Experimental Advanced Superconducting 
Tokamak '°!, except that the gain of the analog signal conditioning component is 8 
V/V to accommodate the signal amplitude of the CLLB detector and PMT output 
(<100 mV), and the FPGA code is different. Data processed by the FPGA algorithm 
are transmitted to the chassis-embedded controller through the PXIe bus in direct 
memory access mode for display or storage °*!. 


The pulse acquisition mode of the FPGA is shown in Fig. 2. L, is the length of 


the entire recorded pulse; L, is the length of the sliding average of the baseline 


calculation; L, is the offset of the pulse trigger, which determines the starting point of 
the baseline calculation; and T, is the moment at which the pulse is triggered. In this 
experiment, the average value of 48—64 ns sampling points before pulse triggering was 


defined as the baseline of the pulse, and the total recording time of the pulse was 1008 


ns. 
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Fig. 2 Pulse recording mode 
2.3 Pulse preprocessing 


First, the detector was energy-calibrated using a '*’Cs gamma source. Then, a 
pulse waveform with an equivalent gamma energy of 3.1-3.3 MeV (the energy band 
containing the thermal neutron peak) was extracted from the **Cf source; therefore, 
we investigated the neutron/gamma discrimination performance in this energy band. A 
total of 33,782 waveforms were measured. Pulse width and amplitude analysis was 
used to remove severely overlapped or chopped waveforms; 33,246 pulse waveforms 
remained after this process. Finally, each pulse baseline value was subtracted from the 


sampled pulse sample data. 


3 Principle of numerical method 
3.1 SK filter 


In 1955, the SK filter based on discrete components was first proposed by R.P. 
Sallen and E.L. Key, and the Gaussian shaping of pulse signals was successfully 
realized. SK filters are widely used for signal filtering and pulse shaping, where they 
outperform other methods in terms of energy resolution and computational workload 
27], The SK filter has a second-order filter circuit with a simple structure and is widely 
used in nuclear signal pulse shaping!**9!, Fig. 3 shows the SK low-pass filter. When 
R1 = R2 =R = 3.99 kQ and Cl = C2 = C = 1.5 pF, the corresponding bandwidth is 
16.7 MHz. 


Fig. 3 Schematic diagram of SK low-pass filter 
3.2 Real-time discrimination method 


The CC method uses the different integration values of the neutron and gamma 
pulse tails as the basis for discrimination. The principle is illustrated in Fig. 4a; the 


discrimination is calculated using Eq. (1): 
Q-Q 
-trL h (1) 
Q Q + Qs 


where Q is the total pulse integral, and Qs and Q, represent the pulse front and tail 


Po 


integrals, respectively. The range of the integration intervals of Qp is determined by ts 
and tm, and the integration interval of Q is determined by t, and te. If the interval 
between t, and t,, is set too large, it will cause a loss of the difference between the 
falling portions of the neutron and gamma pulse signals in Eq. (1), where Pcr thus 
tends to 0. By contrast, if the interval between t, and t,, is set too small, the 


information on useless pulses in Eq. (1) will increase, causing Pc, to tend to 1. 


Therefore, the integration interval needs to be selected according to the pulse 
waveform of the detector response in practical applications. In the CLLB detector, the 
neutron pulse will have smaller Pç; values than the gamma pulse because neutrons 
decay faster than gammas, unlike the typical case of liquid scintillators. 

The principle of the AC method is shown in Fig. 4b, where t, after the peak is 
selected as the moment of amplitude comparison. The different amplitudes of neutrons 
and gammas at this moment are used as the basis for discrimination, where A, and A, 
represent the amplitude of the gamma and neutron pulses, respectively, at the t, 
moment. This method requires the calculation of only one sampling point but is very 
susceptible to the effects of noise. Therefore, the amplitudes were normalized before 


comparison with the other methods. 
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(c) Schematic diagram of time comparison method (d) Schematic diagram of pulse gradient method 


Fig. 4 Schematic diagram of discrimination methods 

The principle of the TC method is shown in Fig. 4c. It can be considered as the 
inverse function of the AC method. A fixed A, threshold line is chosen, and the 
difference in the time it takes for neutron and gamma pulses to decay from their peaks 
to the threshold line is used as the basis for discrimination. T, and T, represent the 
time from the intersection of the decaying portion of the gamma and neutron pulse, 
respectively, with the threshold line A, to the peak. 

The PG method selects peak and post-peak sampling points for gradient 


calculation; the principle is illustrated in Fig. 4d. It can be written as Eq. (2): 


Pg= 
ty — tp 


where Pg is the pulse gradient, A, is the pulse peak amplitude, and A, is the amplitude 


at a specific time interval after the pulse peak. In addition, t, and t, are the 


corresponding moments of the pulse peak and post-peak sampling points, respectively. 


Similarly, the gradient is susceptible to strong amplitude effects and is calculated using 


normalized amplitudes. 


3.3 Evaluation criteria 


The FOM was introduced as an evaluation criterion to objectively assess the 


performance of the discrimination method, as shown in Fig. 5. The FOM is defined °°! 
as shown in Eq. (3): 
D 
FOM = ——___ 
FWHM,„ + FWHM, 
FWHM, and FWHM, are the half-height widths of the neutron and gamma peaks, 


(3) 


respectively, and D is the distance between the neutron and gamma peaks. A larger 


FOM indicates that the method exhibits better neutron/gamma discrimination. 
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Fig. 5 Criteria for evaluating neutron/gamma discrimination 


3.4 MDDM 


Among the four methods described above, only the CC method uses all the 
waveform information. Therefore, it has better noise immunity. The discrimination 


performance of the CC method depends on three parameters, ts, tm, and te. The values 


of t, and t, are usually taken as the moments of the upper edge of the pulse and the 
end of pulse decay, respectively, and the parameter tm is the most difficult to determine. 
We propose a model that maximizes the difference used for discrimination to 
determine the parameter t,, and achieve excellent performance by the CC method. 
When neutrons and gamma rays deposit energy in the CLLB detector, the photon pulse 
generated by the scintillation crystal can be written as shown in Eq. (4). 

I(t) = lye! + Ie" (4) 


where Ty and Ts are the fast and slow components of the decay time constant, 


respectively. Ir and /, represent the proportions of the fast and slow components, 


respectively. 

In the linear region, the pulse output from the CLLB detector can be represented 
by the convolution of the photon pulse with the response function of the PMT and 
readout electronics. However, the final expression of V(t) still contains the factors 
e™t/tf and e—"/*s [31], Because the pulse rises rapidly, this part of the output can be 
replaced by a linear function; thus, the gamma and neutron pulse response described 


above can be written as shown in Eqs. (5) and (6), respectively. 


Kyt ,o< ti 

Vergy Ct) = AÇI pe Etr + Ie ®t"), t > ti (5) 
K,t ‚t< ti 

Veren Ct) = All pne ED + Ipe EN), t> ty (6) 
K, = Kn = A/t, o 


Ify and I,,, represent the fast and slow components, respectively, of the pulse response 


to gamma rays. I, and Isn represent the fast and slow components, respectively, of the 


response to neutrons. A is the response amplitude, which is proportional to the ray 
energy. In the actual detector impulse response, the pulse first rises and then decays, as 
shown in Fig. 6. 
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Fig. 6 Model of output pulse response of CLLB detector 

The rise time of the CLLB response pulse to the peak is the same for rays of 
different energies. tg, t4 and tz represent the moments at which the pulse starts to rise, 
reaches its peak, and is completely decayed, respectively. Są represents the pulse 
integration area of the neutron and gamma pulses from tp to t1. Neutrons and gammas 
of the same energy in the rising phase of the pulse can be considered to have the same 
integration area. S, represents the integration area of the neutron pulse from t4 to tz. 
Because the slow component of the gamma response pulse is more significant in the 
decaying part of the pulse, S, represents the excess integration area of the gamma pulse 
from t4 to tz compared to the neutron pulse. The Pc; values of gammas and neutrons in 


the CC method are calculated as shown in Eqs. (8) and (9), respectively: 


Spr Se — Sys 
P= 8 
OU Si tSp tS. 19) 
Sp — Sng 
Popo — = 


The difference in Pc; between neutrons and gammas is strongly correlated with D 
in Eq. (3). Therefore, the FOM of the CC method is expected to be highest when ô = 
ti +Tm reaches the moment at which the difference in Pey between neutrons and 


gammas is maximum. 


Sns (Sa + Sp + Se) — Sys(Sq + Sp) + SaS 
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(Sa + Sp + Se) (Sa + Sp) 


If Iy < In, Isy > Isn and O< Tm < tz2— t1, Tm is the moment at which the 


(10) 


difference between the Pcr values of gammas and neutrons is largest, as shown in Eq. 


10 


(11). 


TFT I —I 
In = fts n( fnP fya 


11 
Ts— Tf Iq z Isnp ) ( ) 
p = Sa tSp tSoq = Sg tSp p <q (12) 
By calculating T,,, the parameters of the CC method that will yield the best FOM 


can be obtained. 
4 Experimental results and discussion 
4.1 Comparison of methods 


The discrimination effects of the time-domain methods discussed above are all 
related to the parameter values. To objectively evaluate each method, the definition 
domain of the discrimination parameters of each method is first given. Then the best 
FOM using these parameters is obtained as the discrimination performance of the 
method by a trial-and-error approach. Let t, be 10% of the time required for the pulse 
to rise to its peak, t, be the moment of the pulse peak, and t, be 700 ns after t,. For the 
CC method, t,, satisfies Eq. (13); for the AG method, t, satisfies Eq. (14); for the PG 
method, t, satisfies Eq. (15); and for the TC method, the A, range is taken from 0.1 to 
1 peak, at intervals of 1% of the peak. 


t, < tm < te (13) 
A EA (14) 
T (15) 


The FOM of the methods without filtering is shown in Fig. 7. Only the CC 
method can distinguish between neutron and gamma pulses at the current noise level, 
where the parameter tm is taken as 212 ns after t, to obtain the best FOM for this 
method. 
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Fig. 7 Discrimination performance of different methods without filtering 

As show in Fig. 3, the shaping effect of the SK low-pass filter is affected by the 
value of RC. With increasing RC, the cutoff frequency decreases, the attenuation of the 
stop-band increases, and the filtering effect is improved. 

The AC method is more strongly affected by noise, which is strongly correlated 
with the filtering effect. The overall FOM of the TC method is low; it is challenging to 
distinguish neutron and gamma signals correctly because the shaped pulse has a large 
oscillation with time. Therefore, the TC method is not suitable for neutron/gamma 
discrimination. 

The PG method is also more strongly affected by noise. The FOM of the PG 
method is strongly correlated with the filtering effect. In addition, the FOM of the PG 
method is similar to that of the AC method because of the amplitude difference in Eq. 
(3) in the calculation, which indicates that the amplitude difference is the dominant 
factor in discrimination by the PG method. 

Therefore, the CC method is the most suitable processing method for real-time 


neutron/gamma discrimination. 
4.2 Effect of MDDM 


When the t, and tẹ parameters of the CC method are fixed as described in Sect. 
4.1, we take tm = {t|t = 60 + 20+ a,a = 0,1,2,3,...,28} . The FOM of the CC 
method for these values of tf; is shown in Fig. 8. The FOM curve is convex, which 


indicates the best FOM can be obtained by considering the FOM as a function of the 


12 


single parameter t,,, and confirms the suitability of the MDDM proposed in Sect. 3.4. 
Let eg-e; be the Rp region, e;-e2 be the R4 region, and e2-e3 be the Rz region. The 
slope of the curves in the Rọ and R, regions is larger, indicating that the FOM 
improves more rapidly when tm is close to e1. By contrast, the FOM decreases rapidly 
when tm is close to e3. The slope of the curve in the Rz region is smaller, indicating a 
slower decrease in the FOM in this region, although the range is large. Thus, the 
performance of the CC method is somewhat independent of the tm value. 
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Fig. 8 Effect of tm on FOM of CC method 

To solve Eq. (11) for Tm, only If, Isn, Ify, Isy» Tf, Ts and t; need to be calculated. 
They were obtained by fitting the neutron/gamma pulse waveform by Eqs. (5) and (6), 
as follows: 

1. The neutron/gamma pulse signals were extracted separately by the CC method 
(using the parameters that were used to achieve the performance shown in Fig. 7), 
where 2000 pulse waveforms each were extracted for neutrons and gammas. 

2. The average waveforms of neutron and gamma pulses were calculated and 
fitted with Eqs. (5) and (6) for neutrons and gammas, respectively, to obtain 
Tfn TSn, Tfy, and Ts}. 

3. The average of the rise times of the fitted neutron and gamma waveforms is 
obtained as t4 = 30 ns. 

4. The average of Tr, and Ts, is taken as Tf = 81.51 ns, and the average of Tsn 
and Tsy is taken as Ts = 451.27 ns. 

5. By substituting Tp and tT, into Eqs. (5) and (6), we obtained I, = 0.868, Isn = 
0.132, If, = 0.837, and I = 0.163. 
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Fig. 9 Comparison of FOM obtained using MDDM and optimized parameters. The 
black dots represent the potentially optimal FOM for different values of tẹ. The red 
dots represent the FOM obtained using the t,, values calculated using the MDDM for 
different values of te. 


We confirm that the above parameters satisfy Eq. (11). For a fixed value of t,, we 
take te = {t|t = 300 + 20 * a, a = 0,1,2,3, ...,30}. Fig. 9 compares the FOM obtained 
using the t,, values calculated using the MDDM with the potentially optimal FOM for 
different te values. The potentially optimal FOM versus tm curve (black dots) is still 
convex, indicating that te does not need to cover the entire pulse to obtain the best 
discrimination performance of the CC method, and thus the calculation time required 
for integration can be reduced. Experimental tests show that tẹ can be set between 640 
and 740 ns after ts, where the potentially optimal FOM exceeds 1.46. In addition, 
when t, is set in this range, the difference between the FOM obtained using the tm 
value calculated using the MDDM and the potentially optimal FOM is less than 3.9%, 


confirming that it can improve the performance of the CC method. 


5 Conclusion 


This work studied the use of four real-time discrimination methods with CLLB 
detectors. The FOM of the four methods was improved when the pulses were 
adequately filtered. The CC method, with an excellent FOM and noise immunity, is the 
most suitable real-time discrimination method for CLLB detectors. Experimental tests 
showed that when t, in the CC method is set to the moment at which the pulse rises to 
10% of its peak, tẹ should be within 640-740 ns after t, to yield the potentially optimal 
FOM. In this parameter range, the FOM obtained using the tm value calculated using 


14 


the MDDM and the potentially optimal FOM differ by less than 3.9%. This result 
provides a good guide for parameter setting in the CC method. The concept of MDDM 
analysis can be used to model the parameter settings in other time-domain methods and 


provide guidance for parameter selection. 
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